Dynamics and modulation of ring dark soliton in 2D Bose-Einstein 
condensates with tunable interaction 
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We investigate the dynamics and modulation of ring dark soliton in 2D Bose-Einstein condensates 
with tunable interaction both analytically and numerically. The analytic solutions of ring dark 
soliton are derived by using a new transformation method. For shallow ring dark soliton, it is stable 
when the ring is slightly distorted, while for large deformation of the ring, vortex pairs appear and 
they demonstrate novel dynamical behaviors: the vortex pairs will transform into dark lumplike 
solitons and revert to ring dark soliton periodically. Moreover, our results show that the dynamical 
evolution of the ring dark soliton can be dramatically affected by Feshbach resonance, and the 
lifetime of the ring dark soliton can be largely extended which offers a useful method for observing 
the ring dark soliton in future experiments. 

PACS numbers: 03.75.Lm, 05.45.Yv, 03.65.Ge 



I. INTRODUCTION 

Solitons are fundamental excitations of nonlinear me- 
dia and have attracted great interests from diverse con- 
texts of science and engineering, such as dynamics of 
waves in shallow water, transport along DNA and other 
macromolcculcs, and fiber optic communications. The 
realization of Bose-Einstein condensates (BECs) [l| in- 
troduces an unparalleled platform for the study of these 
nonlinear excitations, where both bright and dark soli- 
tons d, 3, H, [1, have been observed. Very recently, 
quasi-lD dark solitons with long lifetime up to 2.8s have 
been created in the laboratory [8| . This offers an unusual 
opportunity to study dark solitons and the relevant the- 
ories. 

Dark solitons are robust localized defects in repulsive 
BECs, which are characterized by a notch in the con- 
densate density and a phase jump across the center 
So far most of the works on dark soliton are limited 
to one-dimensional (ID) BECs ([l3l, for a review see 
[ni, m, [l^l), where they are both stable and easily con- 
trolled in experiments. However, in two dimensions (2D), 
most types of dark solitons are short lived due to dynam- 
ical instability arising from higher dimensionality [j] . For 
example, in ref. [ij] it is shown that both 2D dark lump- 
like soliton and dark stripe soliton will decay into vortex 
pairs under small transverse perturbations. Also these 
2D defects are strongly affected by inhomogeneity of the 
system and suffer from the snaking instability just as in 
optical systems ([§, and references therein). There- 
fore their long-time dynamical behaviors are hard to ob- 
serve in real experiments. 

One candidate for observing the long-time behavior of 
2D dark soliton is ring dark solitons (RDS) which was 
first introduced in the context of optics [H, [l^ and then 



studied in BECs [l8[. In a 2D nonlinear homogeneous 
system, it was first predicted that the instability band of 
a dark stripe soliton can be characterized by a maximum 
perturbation wavcnumber Qmax [l9j ; and if the length of 
the stripe is smaller than the inverse of the wavenumber 
L < 2TT/Qmax, thcn the stripe can be bent into an an- 
nulus with the instability being largely suppressed. This 
was further confirmed for BECs even in an inhomoge- 
neous trap, where both oscillatory and stationary ring 
dark solitons can exist [l^- Besides, the symmetry of 
the ring soliton determines that it is little affected by the 
inhomogeneity of BEC system. Moreover, the study indi- 
cated that the collisions between the RDSs are quasielas- 
tic and their shapes will not be distorted [10]. Due to 
these specific characteristics of the RDS, it has stimu- 
lated great interests on observing 2D dark solitons in 
BECs lUil,!!!,!!!]. Recently, Yang et al. suggested a 
proposal on how to create the RDSs in experiment [2^. 
Nevertheless, the generation of the RDS and their dy- 
namical behaviors in BECs have not been observed in 
real experiment yet. This is because the lifetime of deep 
RDS is not long enough for the experimental observation 
and the stability analysis for shallow RDS has not been 
explored thoroughly. 

In this paper, we first develop a general effective an- 
alytical method to derive the solutions of RDS. This 
is realized by transforming the Gross-Pitaevskii (CP) 
equation with trap potential to a standard nonlinear 
Schrodingcr (NLS) equation. Then we numerically solve 
the 2D CP equation and obtain a stability diagram where 
a stable domain of shallow RDSs can exist. In the un- 
stable region, we find a novel transformation process be- 
tween various dark solitons and vortex pairs. Finally, we 
show that the lifetime of RDSs can be largely extended 
by Feshbach resonance, which is of particular importance 
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for the experimental observation of RDSs. 



II. THE MODEL 

A BEC trapped in an external potential is described 
by a macroscopic wave function ^'(r,<) obeying the GP 
equation (26j , which reads 



ih 



91'(r,i) 
'dt 



where the wave function is normalized by the particle 
number N = J dr \ "if \^ and go{t) = ATTh'^a{t)/m rep- 
resents the strength of interatomic interaction character- 
ized by the s -wave scattering length a{t), which can be 
tuned by Feshbach resonance. The trapping potential 
is assumed to be V{r,z) = m{uj^r^ + w^z^)/2, where 
= .T^ + y^, m is the atom mass, and w^^z are the con- 
finement frequencies in the radial and axial directions 
respectively. Further assuming Q = cUr/tUz <C 1 such 
that the motion of atoms in the z direction is essentially 
frozen to the ground state {f{z)) of the axial harmonic 
trapping potential, the system can be regarded as quasi- 
2D. Then we can separate the degrees of freedom of the 
wave function as ^'(r,t) = tp{x,y,t)f{z), obtaining the 
2D GP equation: 

dib d'^ d'^ Tfi 



where 



Idzlfjz) \' 
Jdzlfiz) P 



It is convenient to introduce the scales characterizing 
the trapping potential: the length, time, and wave func- 
tion are scaled as 

~ , ^ I 
X = ahX,t = — ,^/; = 

LUz ah\/4:Traor] 



respectively, with Uh = y^h/muiz and ao is a constant 
length we choose to measure the time-dependent s -wave 
scattering length. Then the 2D GP equation is reduced 
to a dimensionless form as 

^^ = -iv2^ + 5(^)I^P^ + i^^V^, (1) 

where = d'^/dx^ + d'^/dy-^ = d'^/dr^ + 1/r x d/dr + 
d'^/dO^, il = ojr/oJz, g{t) = a{t)/ao and the tilde is omit- 
ted for simplicity. This is the basic equation we treat 
analytically and numerically. 



III. TRANSFORMATION METHOD AND 
ANALYTIC SOLUTION 

In order to study the dynamics of ring dark soliton, we 
consider the solution of Eq. ([T|) with circular symmetry, 



tp{r,t). In the case of g{t) = C, where C is a nonzero 
positive constant, the main difficulty to solve Eq. ^ is 
the existence of the last term, i.e. the trapping potential. 
Without trap, i.e. = 0, the system is described by the 
standard NLS equation: 



.dQ{R,T) 
dT 



1 d^Q{R,T) 
2^ 9i?2 
C\Q{R,T) 



1 dQ{R,T) 
R OR 
^Q{R,T) = 0, 



(2) 



Under small-amplitude approximation, Eq. ([2]) has been 
transformed to the famous cylindrical KdV (cKdV) equa- 
tion by using the perturbation method fwi . [23j . The 
cKdV equation is known to be basic nonlinear equa- 
tion describing cylindrical and spherical pulse solitons 
in plasmas, electric lattices and fluids (see, e.g., 12711 fo r 
a review) and its exact solution has been derived [2g|. 
Therefore the soliton solutions of Eq. ^ are gained. 

However, there is no effective method to solve Eq. 
(fT|) generally, especially when the last term is time- 
dependent, V(t). Now we develop a method which can 
transform the general form of Eq. ([T]) to Eq. ^ by using 
a transformation: 



V'(r,t) = Q(i?(r,<),r(<))e'''('-'*)+^(*), 



(3) 



where R{r,t), T{t), a{r,t), and c{t) are assumed to be 
real functions and the transformation parameters read 

i?(r,i) a{t)r, 

T{t) = / a^it')dt' + Co, 



(4) 



1 , 



a(r, t) 



1 da{t) 2 
2a{t) dt ^ ' 



where Co is a constant. The condition that such trans- 
formations exist is 



1 d'^a{t) 2 .dait) 2 



Y{t) dt^ a(i)2 dt 



(5) 



Under this condition, all solutions of Eq. ^ can be 
recast into the corresponding solutions of Eq. So 
we build a bridge between the extensive RDS study in 
nonlinear optics (homogeneous system) and the trapped 
BEC system. Furthermore, it is worth while to note that 
the transformation method can be used to solve the gen- 
eral equation: 

du{r,t) ^ ^,^^^ d'^u{r,t) ^ l du{r,t) ^ 
dt ^ ^ '^ dr^ r dr ' 
+ g{t) I u{r, t) p u(r, t) + V{tyu{r, t) = 0, (6) 

when g{t) is proportional to D{t). Eq. ([6]) completely 
describes the dynamics and modulation of both electric 
field in optical systems and macroscopic order parameter 
in atomic BECs in quasi-2D with circular symmetry. 
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In order to get the transformation condition explicitly, 
we substitute y{t) = da/{adt) into Eq. ([5]) and obtain 



This is a standard Riccati equation. We can solve it 
not only for il = constant, but also for various type of 
fi^(t), such as a + bsin{Xt), a + bcosh{t), ae"^* (a, b, A are 
arbitrary constants) and so on. Thus we can study the 
dynamics of system with time-dependent external trap. 
Particularly, when Q is independent on t, the solution of 
Eq. dni) is 

a{t) = Ciscc{nt + C2), (8) 

where Ci and C2 are the integral constants. 

Since under small-amplitude approximation, the solu- 
tions of Eq. ^ have been given, then combining the 
solutions of Eqs. ([2]) and we get the corresponding 
solutions of RDS in the BEC in the external trap poten- 
tial. They are the exact solutions of system when the 
depth of RDS is infinitely small, which help us to under- 
stand the dynamics of RDS. However, when the RDSs 
get deeper, the small-amplitude approximation is invalid 
and we have to appeal to the numerical simulation. In 
the following sections, we study the dynamics and stabil- 
ity of RDS numerically. 

IV. STABILITY OF SHALLOW RING DARK 
SOLITON 

It has been known that starting from the initial config- 
uration with strict circular symmetry, the RDS will os- 
cillate up to a certain time till instabilities develop: shal- 
low RDS slowly decays into radiation and for deep one, 
snaking sets in, leading to formation of vortex-antivortex 
pairs arranged in a robust ringshaped array (vortex clus- 
ter), because of transverse perturbations [I^l- But its 
stability to the perturbation in the radial direction is not 
analyzed. Here we study the stability of shallow RDS 
against the small distortion of the ringshape numerically, 
which is ineluctable in the process of the practical exper- 
iment. 

We study the stability of RDSs by solving Eq. ([T|) 
numerically with the parameter: g{t) = 1, = 0.028 
and the initial radius of RDS Rq — 28.9. Because of large 
initial radius, a reasonable and good approximation is 

V'(x,y,0) = {l-n\y4) 

X [cos^(0)tanhZ(ri) +isin(/)(0)] , (9) 

where r = ^ + y"^, Z{ri) = (ri — i?o) cos (/)(0), ri — 
^(1 — e^)a;^ -I- y^, is the eccentricity of the ring and 
cos (^(0) is proportional to the depth of the input soliton. 
When Cc ^ 0, i?o represents the length of semiminor 
axis of the elliptical configuration. We take Eq. ^ as 
the initial configuration, the validity of which has been 
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FIG. 1: A part of stability diagram of shallow ring dark soli- 
ton with Ro = 28.9 and the trap frequency = 0.028. The 
maximum eccentricity Cmax of the ring in the initial configu- 
ration decreases as the initial depth cos 0(0) of the ring dark 
soliton increases. 



explained in detail and checked in the previous investi- 
gations 

We propagate the 2D time-dependent GP equation 
(Eq. Il])) using two distinct techniques: alternating- 
direction implicit method [2^ . [soj and time-splitting 
Fourier spectral method [sj . The results from these two 
methods are crosschecked and When ec = 0, the prop- 
erties of results are in agreement with those in [l^ very 
well. 

To translate the results into units relevant to the ex- 
periment [1,3, we assume a Rb (ap = 5.7 nm) conden- 
sate of radius 30 fim, containing 20000 atoms in a disk- 
shaped trap with lo^ = 2-11 x 18 Hz and cj^ = 27r x 628 Hz. 
In this case, the RDS considered above has the radius 
i?o = 12.38 /iTO, and the unit of time is 0.25 ms. 

Our numerical results show that the shallow RDS 
(refers to cos</)(0) < 0.67, where the soliton without 
distortion will not suffer from the snaking instability) 
is stable against the small distortion of the ringshape. 
There is a maximum eccentricity emax for a given ini- 
tial depth. When Cc < Cmax, the RDSs are stable and 
oscillate reserving their shape with the same period as 
the unperturbed RDSs until decaying out. The Cmax be- 
comes smaller with the initial depth increasing (see Fig. 
1). This is because lir/Qmax decreases as the RDS be- 
comes deeper [l^, thus the shallower RDS can support 
larger distortion. 

When Cc exceeds Cmaxi snaking sets in and the dark 
soliton breaks into two vortex pairs, presenting a striking 
contrast to the multiples of four pairs reported previously. 
And the evolution of vortex pair is very different from the 
case of the deep RDSs without distortion. 

To illustrate the generic scenarios, we take the typical 
case with cos 0(0) = 0.6 and Cc = 0.4. It initially shrinks 
and when reaching the minimum radius in the short-axis 
direction, it starts snaking and forming two dark lump- 
like solitons in the horizontal direction; they move in the 
opposite direction and then break into two vortex pairs 
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FIG. 2: (Color online) Evolution of the ring dark soliton with 
initial depth cos(^(0) = 0.6, eccentricity Be = 0.4, interaction 
strength g{t) — 1 and the trap frequency Q = 0.028. (a) 
The initial profile shown by color-scale density plots, where 
ring dark soliton corresponds a light-color plot ellipse with 
the length of semiminor axis _Ro ~ 28.9. (b)-(d) correspond 
to i = 80, 400, 540 respectively. The ring dark soliton with 
distortion firstly breaks into two lump solitons, and then two 
pairs of vortex pairs. See text for details. 



(see Fig. 2(b)). The vortex pairs arrange themselves in a 
ring configuration which performs slow radial oscillations. 
Simultaneously the vortices and antivorticcs move along 
the ring (see Fig. 2(c, d)). The result of motion is their 
collision in pairs in the vertical direction, followed by 
vortices and antivortices's mergence and totally forming 
a pair of lumplike solitons (see Fig. 3(a)), which process 
has been confirmed by phase distribution of the system. 
The lumplike solitons move towards to each other, i.e. to 
the center of condensate, in the vertical direction. When 
they reach the minimum radius, a ring dark soliton forms 
(see Fig. 3(b)). After a short time the system returns 
to the state of two dark lumplike solitons. The lumplike 
solitons leave each other (sec Fig. 3(c)) and in the path, 
each of lumplikc solitons breaks into a vortex pair again. 
But the configuration is different to the initial one: the 
vortex (antivortex) is substituted by antivortex (vortex), 
just like the vortex and the antivortex passing each other 
directly following the motion before the mergence. Then 
the vortices and the antivortices keep on moving (see Fig. 
3(d)), then collide and merge in the horizontal direction. 
The dynamical process repeats itself. 

This dynamical behavior of the vortex pairs is novel 
and part of the process is similar to the behavior of lump 
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FIG. 3: (Color online) Evolution of vortex pairs following 
the Fig. 2 for (a) t = 640, (b) t = 660, (c) t = 680 and 
(d) t = 1000. The vortex pairs collide and merge into dark 
lumplike solitons, which move, then form a ring dark soliton 
for a short time and revive. Finally the lumplike solitons 
again break into vortex pairs. See text for details. 

soliton reported in [TJ] ■ The certain centrical collision of 
vortex pair provides a potential tool to study the collision 
dynamics of vortices. 

V. CONTROLLING OF DEEP RING DARK 
SOLITONS 

The deep ring dark solitons (refer to cos(/)(0) > 0.67) 
suffer from the snaking instability and can only survive 
for a short time (typically 10 ms in the numerical simu- 
lation) before changing into other soliton type. Longer 
lifetime is necessary for the practical observation of com- 
plex soliton physics such as oscillations or collisions [1]. 
Our results show that the Feshbach resonance manage- 
ment technique can largely extend the lifetime of the ring 
dark solitons, which makes it possible to study the long- 
time behavior of 2D dark soliton in experiment. 

Feshbach resonance (3^ is a quite effective mechanism 
that can be used to manipulate the interatomic inter- 
action (i.e. the magnitude and sign of the scattering 
length), which has been used in many important exper- 
imental investigations, such as the formation of bright 
solitons 0. Feshbach resonance management refers to 
the time-periodic changes of magnitude and/ or sign of 
scattering length by Feshbach resonance [3^, and has 
been widely used to modulate and stabilize bright soli- 
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FIG. 4: (Color online) The snapshots of condensate with ring 
dark soliton under the Feshbach resonance management. The 
parameter are given as follow: (a) Ec = 0, cos (^(0) = 1, g{t) — 
2 - sm{2nt) and f = 60 (« 15 ms); (b) Cc = 0, cos 0(0) = 1, 
g{t) — e^'^* and t — 100 (~ 25 ms). Other parameters are 
the same as Fig. 2. In the former case the ring dark soliton 
will break into 20 vortex pairs, while only 4 vortex pairs in 
the later one. The total atom number is the same in these 
two plot, so for clearness, we use different color scale. 



When the system subjects to the modulation of g{t) = 
1 — sin(ilt), the RDS can exist for longer time before 
breakup. For example, when cos(^(0) = 0.76, the life- 
time can be extended up to 45 ms from 10 ms, in which 
process, we can observe one complete cycle of oscillation 
of RDS. While if we just reduce the scattering length but 
keeping it constant, the lifetime of RDS will not change 
much. For example, if g{t) = 0.2, the RDS will become 
snaking at about 15 ms. So the lifetime extension effect 
is due to the Feshbach resonance management. From Eq. 
([5]), we can see that fl is one intrinsic frequency of sys- 
tem. So it is reasonable to deduce that this effect is a 
resonance phenomenon. Furthermore, this effect is even 
valid for shallow RDS suffering from the instability due to 
the distortion of ringshape. It has been check that under 
Feshbach resonance management g{t) = 1 — sin(f7t), the 
RDS with cos (f>{0) = 0.6 and ec = 0.4, has a lifetime of 50 
ms, which is much larger than the original 10 ms. Thus 
under large perturbation the RDS can exist long enough 
to be observed, which make the experimental study of 
RDS easier. 



tons [3J,l35|, while its effect on the dynamics and stability 
of dark soliton is not very definite [36| . Taking the ring 
dark soliton as an example, we demonstrate that the Fes- 
hbach resonance can dramatically affect the dynamics of 
dark solitons. 

To show the effect of Feshbach resonance management, 
Eq. ([T|) is integrated numerically with eccentricity Cc = 
and a time dependent nonlinear coefRcient g{t); other 
parameters and the initial condition are the same as in 
Sec. IV. We consider the cases of git) - ae^'^* and 
g{t) ~ a ± sin(w<), where a and uj are arbitrary con- 
stant parameters. Our results show that for the deep 
ring dark solitons, the Feshbach resonance management 
remarkably changes the evolution and the instability of 
solitons, indicated by the number of vortex pairs which 
arise due to the snaking instability. Two typical exam- 
ples are shown in Fig. 4: in the former one, the black 
ring dark solitons (refer to cos 0(0) = 1) break into 20 
vortex pairs, while the latter one only has 4 pairs. These 
phenomena are very different from those with constant 
scattering length which have 16 vortex pairs. This makes 
it possible to study the dynamics of snaking instability 
in detail. What is more interesting and more important 
are the cases in which the lifetime of the RDS can be 
extended, such as g{t) = 1 — sin(rit), as discussed below. 



VI. CONCLUSION 

In conclusion, we have introduced a new transforma- 
tion method to derive the solution of ring dark soliton, 
which provides a powerful analytic tool to study the 2D 
BECs and nonlinear optical systems with circular sym- 
metry. Then we check the stability of shallow ring dark 
soliton, and find a novel dynamical behavior mode of 
solitons. Furthermore, we study the effect of Feshbach 
resonance management on the evolution and the stabil- 
ity of ring dark soliton, also discover a method to extend 
the lifetime of ring dark soliton largely. We show that the 
ring dark solition is a potential candidate for observing 
the long-time behavior of 2D dark soliton. 
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